Hole spin mixing in InAs Quantum Dot Molecules 
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Holes confined in single InAs quantum dots have recently emerged as a promising system for 
the storage or manipulation of quantum information. These holes are often assumed to have only 
heavy-hole character and further assumed to have no mixing between orthogonal heavy hole spin 
projections (in the absence of a transverse magnetic field). The same assumption has been applied 
to InAs quantum dot molecules formed by two stacked InAs quantum dots that are coupled by 
coherent tunneling of the hole between the two dots. We present experimental evidence of the 
existence of a hole spin mixing term obtained with magneto-photoluminescence spectroscopy on 
such InAs quantum dot molecules. We use a Luttinger spinor model to explain the physical origin 
of this hole spin mixing term: misalignment of the dots along the stacking direction breaks the 
angular symmetry and allows mixing through the light-hole component of the spinor. We discuss 
how this novel spin mixing mechanism may offer new spin manipulation opportunities that are 
unique to holes. 

PACS numbers: 78.20.Ls, 78.47.-p, 78.55.Cr, 78.67.Hc 



I. INTRODUCTION 

An electron excited across the band gap of a semicon- 
ductor quantum dot (QD) leaves behind a hole in the 
otherwise full valence states. This hole behaves like a 
charged particle, much like the electron, though with a 
substantially larger effective mass. Because the valence 
states are derived from p-type atomic states of the lat- 
tice, a hole experiences a strong spin-orbit (SO) inter- 
action that leads to a new spin basis in which the low 
energy hole state has total angular momentum J = 3/2. 
Because light hole (LH) states {Jz ~ ±1/2) are shifted to 
larger energy by confinement and strain, it is generally a 
good approximation to treat low-energy holes in QDs as 
if they have only heavy-hole (HH) ( = ±3/2) character 
with a pseudo-spin ^ . This simple picture has proven to 
be remarkably useful and explains a wide variety of op- 
tical and magneto-optical properties of quantum dotSfii^ 
In the absence of a transverse magnetic field, the orthog- 
onal HH spin projections do not mix. The well isolated 
spin projections and lack of a significant contact hyper- 
fine interaction with nuclear spins makes the HH spin 
projection a good candidate for the storage of quantum 
information.i^ 

Here we present experimental evidence of mixing be- 
tween HH spin projections in coupled quantum dots. The 
mixing is observed in the optical spectra of stacked pairs 
of self-assembled InAs quantum dots near the applied 
electric field that induces coherent hole tunnelling be- 
tween the QDs. Spin mixing between bright and dark 
exciton spin configurations causes dark states to gain op- 
tical intensity, which we have observed in a number of 
cases. Here we present an example in which the spin 
mixing is sufficiently large that we can directly observe 
anticrossings between bright and dark exciton states and 
measure the mixing between opposite heavy-hole spin 



projections. In Sect. [IT] we present the energy levels for 
the neutral exciton states of a coupled pair of QDs and 
describe the expected bright and dark state energy lev- 
els and interactions in the absence of hole spin mixing. 
In Sect, mil we present experimental evidence of the ap- 
pearance of bright-dark anticrossings and show that the 
experimental data can be phenomenologically explained 
by the presence of hole spin mixing. In Sect. IIVI we use 
a Luttinger spinor model to explain how the spin mix- 
ing can arise when misalignment of the QDs along the 
stacking axis breaks the molecular symmetry. A com- 
plete description of the theoretical model can be found 
in App.|B] In Sect. |V] we summarize our observation and 
explanation of the origin of hole spin mixing. 



To build a quantum information processing de- 
vice around hole spins, it is necessary to have op- 
tical or electrical mechanisms for coherently creating 
and manipulating superpositions of orthogonal spin 
projectionsj^ii The novel spin-mixing mechanism pre- 
sented here presents new opportunities for electrical con- 
trol of hole spin projections, which are discussed in 
Sect. I VII We calculate the purity of the hole spin states 
as a function of applied electric field and show that the 
heavy-hole states remain good states for the storage of 
quantum information away from the electric field of tun- 
nel coupling. We suggest a possible scenario for electri- 
cally gating the hole spin mixing and thus manipulating 
the hole spin projection. We note that an electron spin- 
mixing anticrossing was previously measured and used 
for optical spin controli^ The hole spin anticrossing en- 
ergy measured here is an order of magnitude larger. 
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II. ENERGY LEVEL STRUCTURE OF THE 
NEUTRAL EXCITON 

When two InAs QDs are stacked on top of one another, 
electrons or holes can tunnel between the two dots to 
create quantum dot molecules (QDMs)i^ In general the 
tunnel coupling is weak because the natural distribution 
of dot size, shape, and alloying leads to different confined 
energy levels in each dot. An electric field applied along 
the growth direction can tunc the energy levels into res- 
onance to enable coherent tunneling of electrons or holes 
between the dots and the formation of delocalized states 
with molecular orbital character The experimental sig- 
nature of the delocalized molecular orbitals is the forma- 
tion of an anticrossing between photoluminesccncc (PL) 
lines that come from direct (electron and hole in the same 
quantum dot) and indirect (electron and hole in differ- 
ent quantum dots) states. Fig. [TJi shows a calculation of 
the energy levels and anticrossings for the neutral cxci- 
ton (X" : one electron and one hole) in a QDM where the 
dots arc separated by a 4 nm barrier and the hole tunnels 
between dots. Direct and indirect states are labeled. 
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FIG. 1: (Color Online) a) Energy levels of bright and dark 
exciton states in a QDM at B = T. b,c) Energy levels of 
bright states (b) and bright and dark states (c) at B = 6 
T. Black (red) lines indicate optically bright (dark) exciton 
states. 

To denote the spatial locations of the electrons and 
holes we use the notation (^^^^), where cb (gt) are the 
spins of the electrons in the bottom (top) dot and simi- 
larly for holes. This notation describes the sates far away 
from an anticrossing. At the anticrossings, the molecu- 
lar states can be described as symmetric and antisym- 
metric combinations of these basis states. If we make 
the usual assumption that holes have only HH char- 
acter {J. = ±3/2), there are four possible pairings of 
the electron and hole spin projections in a neutral exci- 
ton. We use 1, I to indicate the electron spin projection 
(5e = ±1/2) and ff, 4 to indicate the HH spin projec- 
tion (Jz = ±3/2). The two spin pairings "f-tj- and have 
exciton angular momentum ±1 and can couple to pho- 
tons. These exciton spin configurations are called bright 
states. The other two spin configurations (ttT and |-1|) 
have exciton angular momentum ±2 and are called dark 
excitons because they do not couple to photons. One set 
of bright and dark states are labeled in Fig. [1^ using our 



notation. The direct bright and dark states are split due 
to the isotropic electron-hole exchange interaction. This 
interaction is suppressed when electrons and holes arc in 
separate dots and thus the indirect bright and dark states 
are degenerate. In a typical experimental spectra, only 
the bright states are evident. 

In the absence of a magnetic field or hole spin mixing, 
orthogonal electron and hole spin projections are degen- 
erate. In this work we ignore the small mixing between 
the two bright states that arises from the anisotropic ex- 
change interactionfii The two bright exciton spin config- 
urations are thus degenerate and only the electron spin 
down configuration is labeled in Fig.[T^. One anticrossing 
is observed where the direct and indirect bright exciton 
spin configurations would be expected to intersect (the 
upper anticrossing in Fig. [T^). Another anticrossing is 
observed for the dark exciton states (lower anticrossing). 
Tunnel coupling, anticrossings, and spin interactions of 
electrons and holes in both neutral and charged excitons 
have been observed fi^iHiiiii^ In this work we focus on 
the tunnel coupling of holes in the neutral exciton state, 
as depicted in Fig. [TJi. 

When a longitudinal magnetic field is applied (Faraday 
geometry: parallel to the optical axis), the natural expec- 
tation is that a Zeeman splitting between the two bright 
exciton configurations will be observed, with a splitting 
proportional to the sum of the electron and heavy hole 
g factors. However, in QDMs where holes tunnel the 
formation of molecular orbitals substantially alters the 
hole g factor i^ii^iii Consequently, the hole g factor de- 
pends on the applied electric field and strong resonant 
enhancement or suppression of the Zeeman splitting is 
observed at the electric field of couplingF^^ 

In Fig. [IId we show the calculated energies of the bright 
states from Fig. [Iji when a magnetic field of 6 T is ap- 
plied. The first effect of the magnetic field is a Zeeman 
splitting of the two bright states that were degenerate in 
Fig. [1^. This is seen most clearly at the edges of Fig.[Tl3: 
the degenerate lines at +0.1 meV in Fig. [IJi split and 
move to -0.3 meV and ±0.5 meV. The second effect of the 
magnetic field is the introduction of a g factor resonance. 
The g factor resonance is seen most clearly by looking at 
the anticrossings. The degenerate anticrossings in Fig.[TJi 
have an anticrossing gap of 214 ^eV. In Fig. [T|d, the an- 
ticrossing gap for the lower Zeeman branch expands to 
approximately 400 /icV, while the anticrossing of the up- 
per Zeeman branch collapses to approximately 30 ficV. 

The difference in anticrossing energies arises because 
different molecular orbitals have different g factor con- 
tributions from the barrier, which determine the net 
hole g factorJ^iii The two lower energy molecular states 
(one for each bright electron-hole spin orientation) have 
antisymmetric (noded) orbital character (lower inset to 
Fig. [TJd) . The node in the molecular wavefunction sup- 
presses the contribution of the GaAs barrier to the net 
hole g factor. Because the barrier hole g factor is pos- 
itive, suppression of this contribution increases the rel- 
ative weight of the negative g factor contribution from 
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the InAs QDs. The increase in the magnitude of the g 
factor (more negative) causes the Zeeman spHtting of the 
two antisymmetric (noded) molecular orbitals to increase 
on resonance. Conversely, the higher energy molecular 
states have symmetric (node-less) orbital character (up- 
per inset). The contribution of the barrier is thus en- 
hanced on resonance, offsetting the negative contribution 
from the InAs QDs and reducing the Zeeman splitting of 
these two lines. The combination of the enhanced Zee- 
man splitting for one molecular branch and suppressed 
Zeeman splitting for the other branch leads to the differ- 
ent anticrossing energies. In Fig. [TJ: we plot the energies 
of both the bright and dark exciton states including the 
Zeeman splitting and the resonant changes in g factor. In 
the absence of hole spin mixing there are no anticrossings 
where bright and dark states intersect. 

The counterintuitive antisymmetric (noded) character 
of the molecular ground states in Fig.fTJj is a consequence 
of the spin-orbit interaction, which mixes HH and LH 
statesiii When the barrier separating the QDs is thin, 
the molecular ground state has bonding (symmetric) or- 
bital character and the first molecular excited state has 
antibonding (antisymmetric) character, in analogy with 
natural diatomic molecules. As the thickness of the bar- 
rier is increased (to 4 nm in Fig. [T]), the contribution of 
the LH states becomes more important and leads to the 
reversal of the orbital character. The reversal provides 
one indication that LHs can not be neglected in QDMs. 
To include LHs, hole states are described as Luttinger 
spinors that contain all four projections of Jz, but each 
spinor is dominated by a single HH spin projectioUfi^ As 
we describe below, the mixing that leads to the reversal 
of symmetric and anti-symmetric orbital states does not 
result in mixing between these HH spin projections and 
consequently cannot explain the appearance of bright- 
dark anticrossings. However, if misalignment of the dots 
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So is the electron-hole exchange interaction that splits 
bright and dark states when the electron and hole are 
in the same dot. //_b is the Bohr magneton, B is the 
magnetic field, is the electron g factor in the bottom 
dot. ghB {ghr) is the g factor for a hole in the bot- 
tom (top) dot. txo is the tunneling matrix element, d 
is the effective barrier thickness, which determines the 



along the stacking axis breaks the symmetry of the QDM, 
at electric fields near the point of tunnel coupling the 
spin-orbit interaction combines with the broken symme- 
try to permit mixing of spinors with different HH spin 
components. 

Fig [2^ and d display calculated spectra of the anti- 
crossings of the X'^ state at B = and 6 T. These cal- 
culations use the same parameters as Fig. [U but include 
a color mapping to display the optical intensity of the 
lines. The calculated spectra in Figl^Ji. and d are rep- 
resentative of the typical experimentally observed spec- 
tra. To calculate these spectra we use matrix Hamil- 
tonians, which have been shown to provide an accurate 
phcnomenological model of tunneling, spin interactions, 
and resonant changes to g factor in a wide variety of 
samples] The basis states are the possi- 
ble spatial and spin distributions of the electron and hole. 
When an electric field is applied, the lowest electron en- 
ergy level in the top dot is at significantly higher energy 
than the confined electron energy level of the bottom dot 
and the states with the electron in the top dot can be 
neglected. We can describe the basis states in our no- 
tation {{I^iIt)^ ^^^'^ holes denoted by the dominant HH 
spin projection. Using this notation, the basis states of 
the Hamiltonian are: 



(lo\ fio\ fio\ fio\ 

Note that only the electron spin down case is shown 
because the overall matrix is block diagonal for the two 
electron spin projections. The first two states are bright 
excitons, the second two dark excitons. The Hamiltonian 
that describes the energy of the neutral exciton state is: 



hm 



(2) 



I 

slope of the indirect lines when the electric field, F, is 
applied, hm is the hole mixing term, which is set to zero 
for the calculations in Fig|2^ and d. gi2 is the resonant 
contribution to the g factor from the barrier, which has 
opposite sign to ghr^ The relative signs of txo and gi2 
insure that the bonding orbital has a decreased splitting. 
Here txo < because the molecular ground state has 
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antibonding character. 

The energies of the neutral exciton states are calcu- 
lated by finding the eigenvalues of the matrix at a spe- 
cific value of the field, F. Because the final state after 
optical recombination contains no particles, the energies 
of the neutral exciton initial states are exactly the ener- 
gies of the observed PL lines. The optical intensities are 
calculated by multiplying the corresponding eigenvector 
by an optical intensity vector: 



/ 1 \ 

/nd/^/2 






(3) 



which simply gives unit intensity to the direct transi- 
tion (electron and hole in the same dot) and a fraction 
of that intensity (Ind) to the indirect transition (elec- 
tron and hole in different dots). The dark states have 
no optical intensity. In the calculations, therefore, any 
optical intensity for dark states must come from mixing 
with bright states. 
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FIG. 2: (Color Online) Experimental (b) and calculated (a,c) 
photoluminescence spectral map of the neutral exciton at B 
= T without (a) and with (c) phenomenological hole spin 
mixing term, (d-f) Spectral maps as in (a-c) with B = 6 T. 
Scales for panels (a-c) are the same and are indicated in panel 
a, similarly for (d-f). Callout numbers are referenced in the 
text. 



III. EXPERIMENTAL EVIDENCE OF HOLE 
SPIN MIXING 

The calculated PL spectral maps in Figl2^ and d are 
representative of the observed behavior for the anticross- 
ing of the neutral exciton in most samples. As mentioned 
in Sect. HI we have observed a number of QDMs where the 
dark states gain optical intensity in the vicinity of the an- 
ticrossing region. This intensity gain cannot be explained 
with any of the previously observed QDM properties, in- 
cluding spin-conserving tunneling and molecular g factor 
resonances. Figj^b and e shows experimental data for a 
QDM where the dots are separated by a 4 nm tunnel bar- 
rier. In this example, not only do the dark states gain op- 
tical intensity, we also observe new anticrossings between 
bright and dark states at high magnetic field. These an- 
ticrossings allow us to directly measure the magnitude 
of the spin mixing and to show that it is responsible for 
the dark state intensity. In the zero magnetic field case 
(Figl^b) the new anticrossing with the dark exciton state 
is indicated by callout 1. At B = 6 T (Fig[2};) a complex 
pattern of additional anticrossings (callouts 3-6) appear 
near the electric field of the tunnel resonance. We will 
show that all of the additional anticrossings can be ex- 
plained by the inclusion of the hole spin mixing term, 
which couples bright and dark exciton states. 

The calculations in Fig[2^ and d are obtained with nu- 
merical values determined by the experimental data in 
Figl2b and e with the hole spin mixing term (hm) set to 
zero (see Appendix |A|. It is clear that when the hole spin 
mixing term is set to zero the phenomenological Hamil- 
tonian (Figl2^ and d) does not capture all of the features 
that appear in the experimental data. However, as shown 
in Figllfc and f, all of the dark states and anticrossings 
in the experimental spectra are explained if we turn on 
the hole mixing term. Best agreement with the data is 
obtained when hm = 92/xeV, tx" = and all other pa- 
rameters have the same value as in Figl2^ and d. The 
suppression of tx" =0 is discussed below. 

Without the hole spin mixing term, bright and dark ex- 
citon configurations are independent. In this sample the 
intradot (direct) dark exciton configurations typically lie 
about 200 fieV below the bright exciton configurationsi^ 
The anticrossing observed in Figl2k occurs between the 
direct and indirect bright exciton states. Spin-conserving 
tunneling also couples the direct and indirect dark exci- 
ton states, but no signature appears in the PL spectra 
because the dark states do not couple to photons. The 
hole spin mixing term allows the spin-up hole in the bot- 
tom dot to mix with the spin-down hole in the top dot. 
This coupling mixes bright (e.g. Itf) and dark (e.g. j-lj.) 
exciton states and creates anticrossings wherever bright 
or dark exciton states would cross. The eigenstates with 
mostly dark exciton character gain optical intensity near 
the anticrossings because they contain nonzero bright ex- 
citon components as a result of the mixing. 

In Figl^b callout 1 points out the dark exciton states 
that have gained optical intensity in the experimental 
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spectra. In Fig[2}: callout 2 indicates that the inclusion 
of a the hole spin mixing term in the phenomenological 
Hamiltonian leads to the dark exciton states gaining op- 
tical intensity. The direct dark exciton state ( ) is no 

longer an cigcnstatc of the system. As a result of the hole 
spin mixing, the new eigenstate includes a nonzero con- 
tribution from the indirect bright exciton state ((^o°))- 

It is this bright component that gives the eigenstate op- 
tical intensity. The appearance of the dark state at zero 
magnetic field can not be explained without the hole spin 
mixing term. 

When a magnetic field is applied, both bright and dark 
exciton configurations undergo a Zecman splitting. In 
the absence of hole spin mixing the bright and dark states 
simply cross and the dark states remain dark (Fig[2]i). In 
the presence of hole spin mixing, each crossing of bright 
and dark states becomes an anticrossing observable in 
the experimental spectra. Callout 3 in Fig. [5^, for exam- 
ple, is an anticrossing between a direct dark and indirect 
bright exciton state. We can directly measure the magni- 
tude of this anticrossing gap (180 fxeV) to determine the 
magnitude of the spin mixing term. Fig. [2f shows that 
the inclusion of the spin mixing term explains the addi- 
tional anticrossings observed in the experimental spec- 
tra. The anticrossings indicated by callouts 7 and 8 arise 

from the anticrossing of the indirect bright ((^o^^) a-nd 

direct dark ((^|'o)) cxcitons. These calculated anticross- 
ings correspond to the observed anticrossings indicated 
by callouts 3 and 4 in Fig.[2]3. The anticrossings indicated 
by callouts 9 and 10 arise from anticrossings between the 

direct bright ((^|'o)) ^^^^ indirect dark (^q'5|^) states 
and correspond to the observed anticrossings indicated 
by callouts 5 and 6 in Fig. [2^. The explanation of all of 
these complex anticrossing patterns by the inclusion of a 
single term in the matrix Hamiltonians provides strong 
phenomenological evidence for the existence of hole spin 
mixing. In the next section we will address the physical 
origin of such a mixing term in the misalignment of QDs 
along the stacking axis. 

The phenomenological matrix Hamiltonians we use 
have a limited capacity to make quantitative predictions 
in the case of hole spin mixing. In typical QDMs that 
do not show hole spin mixing, it is possible to measure 
each parameter independently in order to construct a 
quantitatively reasonably matrix Hamiltonian or to fit 
a single undetermined parameter. In this case, the hole 
spin mixing is of the same order of magnitude as spin- 
conserving tunneling, electron-hole exchange and the res- 
onant change in g factor. Consequently it is impossi- 
ble to determine each parameter independently or to ob- 
tain a quantitative fit to all parameters simultaneously. 
This limitation manifests in the suppression of the spin- 
conserving tunneling rate txo ■ Lateral offset between 
the QDs is expected to suppress the tunneling rate and 
the reversal of molecular orbitals predicts that the tun- 



neling rate should be very small for barrier thicknesses 
near those of this sampleJ^ A very small tunneling rate 
is therefore entirely plausible but the limitations of our 
model do not allow us to be more quantitative. 

We conclude this section with a discussion of the ex- 
tent of the experimental evidence for the presence of the 
hole spin mixing term. We have observed PL patterns 
evidencing hole spin mixing in the neutral exciton spec- 
tra of 6 other QDMs for samples with barrier thicknesses 
of 3 and 4nm. Preliminary work indicates that the effect 
is much smaller in a sample with a 6nm barrier. The 
experimental data presented in Figl2j3 and e provides a 
striking and clear example of the contributions of the hole 
spin mixing term because the tunneling term is unusually 
small. This small tunneling rate is consistent with, but 
not proof of, a lateral misalignment between dotsjSi 

We find additional experimental confirmation of the 
existence of a hole spin mixing term in the magneto-PL 
spectra of the positively charged trion (not shown) . The 
excited state of the positive trion contains two hole spins, 
which can be in a triplet configuration when the holes are 
in separate dots. The hole spin mixing again introduces 
new anticrossings in the experimental spectra, in this case 
between the triplet states and singlet states that have 
both holes in the same dot. The additional anticrossings 
that appear in both the neutral exciton and positive trion 
spectra are enabled by the hole spin mixing term, which 
allows spin-flip tunneling. The spin-flip tunneling can 
be seen in the phenomenological Hamiltonian (Eq. 
the hole spin mixing term connects eigenstates that have 
holes in different dots with different spin orientations. 



IV. PHYSICAL ORIGIN OF THE HOLE SPIN 
MIXING TERM 

Although we can model the observed data with the 
addition of a phenomenological spin mixing term, the 
question of the physical origin of such a term remains. 
In analogy with single QDs, one might suspect Rashba 
and Dresselhaus SO interactions.i^i^i However the large 
magnitude of the spin anticrossing gap in Fig. [21 A ~ 180 
/ieV, is not consistent with the high spin purity reported 
for holes in single QDsi^i^ Instead, it seems that the 
QDM geometry has enabled a new mechanism leading to 
strong spin mixing. We propose that such a mechanism is 
the SO interaction mediated by LH. In bulk semiconduc- 
tors LH are known to couple HH states with orthogonal 
spinsi^ This effect is small in single QDs because LH 
are high in energyi^ In QDMs, however, the small ef- 
fective mass of LH causes them to have large tunneling 
rates. As a result, bonding LH states are close in en- 
ergy to the lowest-lying HH statesj^^ and their influence 
becomes important. 

To study the effect of the valence band SO interac- 
tion, we use the simplest description of hole states in- 
cluding HH-LH coupling: the four-band Luttinger-Kohn 
HamiltoniaufSS. The solutions of this Hamiltonian are 



FIG. 3: Schematic depiction of lens-shaped QDs in an ideal 
QDM with no misalignment (a) and a QDM with misalign- 
ment (b). c) Cross-sectional STM image of vertically stacked 
In As QDs showing misalignment. 



Luttinger spinors, four component objects with two HH 
and two LH components. The expression of a Luttinger 
spinor in an ideal QDM, formed by two identical lens- 
shaped QDs perfectly aligned along the stacking axis (left 
panel of Fig. E]) is-J^^ 



/c+3/^Jr)|J, = +§) 



+ 1 ./m,+2(r) \Jz = -\ 
3 
2 



(4) 



Here | > is the Bloch part of the wavefunction, /(r) 
is the envelope function, and cj, a numerical coefficient 
which gives the weight of each component. The enve- 
lope components of the spinor have the symmetries of 
the confining potential. Since the ideal QDM has circu- 
lar symmetry, we can label each of the components by 
an envelope angular momentum = 0, ±1, ±2 . . .. The 
complete Luttinger spinor, however, does not have circu- 
lar symmetry because it is broken by the valence band 
SO interaction. Instead, the spinor can be classified by 
the total angular momentum Fz = mz + Jz and the main 
quantum number k. 

As can be seen in Eq. Q , the spinor contains a mixture 
of tr [Jz = +3/2) and 4| {Jz = -3/2) HH components. 
Despite this mixture, the low-lying hole states of a QD 
(1^2 = ±3/2)) arc dominated by the HH with = 0. 
In typical InAs QDMs, one HH spin component makes 
up over 95% of the ground state weight, with small con- 
tributions from the LH components (less than 5%) and 
the HH component with opposite spin (less than 0.1%). 
One can then identify the \Fz = +3/2) and \Fz = -3/2) 
spinors with the f|- and ^ HH of the usual single-band 
description. 

To show that the weak mixture of HH spins within the 
spinor of a circularly symmetric QDM is not responsible 
for the features observed in Fig. [21 we calculate the low- 
energy hole states of an ideal InAs/GaAs QDM subject to 
a longitudinal magnetic field of S = 6 T (see Appendix iBl 
for calculation details) . The resulting energy spectrum is 
plotted in Fig. 2^. The \Fz,k) symmetry of the states is 
indicated, along with the spin of the dominant HH com- 
ponent of each spinor. In the figure we see anticrossings 



01 

o 

X 



-3/2,2>J 



104 - 




105 




-4-2 2 

F (kV/cm) 



-4-2 2 

F (kV/cm) 



FIG. 4: (Color online). Hole energy levels vs. electric field at 
B = 6 T, calculated with k-p theory: (a) no misalignment; (b) 
misaligned QDM with 3 nm lateral offset. In (a) we indicate 
the \Fz, k) quantum numbers and the pseudo-spin. Note that 
the states with opposite spin cross (anticross) in the absence 
(presence) of misalignment. The inset in (b) shows the mag- 
nitude of the spin anticrossing gap as a function of the lateral 
offset for various interdot barriers. 



at resonant electric fields {Ez ^ kV/cm), where states 
with the same Fz mix to form bonding and antibonding 
molecular states. These anticrossings correspond to the 
spin-conserving tunneling observed in typical experimen- 
tal spectra. In addition. Fig. shows crossings between 
levels with different Fz (different pseudo-spin). These 
states cross (rather than anticross) because the SO inter- 
action does not mix states with different Fz . These cross- 
ings, however, occur between the states that anticross in 
both the experimental spectra and the phenomenologi- 
cal calculations that include the spin mixing term hm 
(Fig. [2]). The absence of anticrossings at these points in 
Fig. 2^ demonstrates that the inclusion of the SO inter- 
action is not sufficient to explain the new experimentally 
observed anticrossings. 

In order to obtain anticrossings that match the ex- 
perimental spectra, the Fz symmetry must be broken. 
The total angular momentum symmetry can be removed 
by structural distortions breaking the circular symme- 
try of the QDM. The distortion could be dot eccentric- 
ity, which is often present in Stranski-Krastanov grown 
QDsi^ However, our simulations (not shown) indicate 
that eccentricity only weakly mixes states with different 
Fz , and the anticrossings gaps it produces never reach the 
large experimental value. Instead we consider a lateral 
offset between the QDs which form the QDM, as shown 
in Fig. ^p. QDM misalignment efficiently removes the 
circular symmetry at resonant electric fields, leading to a 
strong mixing of states with different Fz and, as we show 
below, to spin anticrossing gaps comparable to those of 
the experiment. 

The effect of misalignment is illustrated in Fig. 
where we plot the hole energy spectrum for the same 
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QDM as in Fig. [4^, but now including a lateral offset 
of 3 nm. One can see that the presence of misalignment 
introduces the expected anticrossings between states with 
opposite pseudo-spin. 

The inset in Fig. |4|d shows the magnitude of the an- 
ticrossing gap between states with different pseudo-spin, 
A, as a function of the lateral offset. InAs QDMs with 
different interdot barrier thickness are considered.'''^ A 
increases linearly with the offset, and for 5 nm it may 
reach values of 200 /ieV, which are comparable to the 
value observed in the magneto-photoluminescence spec- 
tra of Fig. O The inset also reveals that the effect of 
the barrier thickness (d) and height (valence band off- 
set) is important. In general, the weaker the tunneling 
the smaller the anticrossing gap (compare e.g. d=1.7 nm 
and d=6 nm). This is because the bonding LH states 
are farther in energy, and their influence decreases. The 
nature of the hole molecular state is also relevant. At 
d=1.7 nm the calculated ground state is bonding, but 
it switches to antibonding at d=1.8 nm. The antibond- 
ing ground state contains a larger admixture of spinor 
components'^, which explains the drastic increase in an- 
ticrossing gap as a result of a small increase in barrier 
thickness. 

To further support misalignment as the origin of the 
phenomenological spin mixing term of Eq. ^ we show 
that the inclusion of lateral offset in the k-p theory in- 
troduces anticrossings between exactly the same states 
that anticross when the spin mixing term is included in 
the phenomenological Hamiltonian. In Fig. [5] we com- 
pare the exciton emission spectrum calculated with the 
phenomenological Hamiltonian (top row) and k-p theory 
(bottom row). The parameters of the phenomenological 
Hamiltonian are the same as in Fig. [21 but the resonant g- 
factor and electron-hole exchange have been neglected for 
simplicity. Because these terms have been neglected, the 
results should not be compared with Fig. [Jf. Fig. [51(a) 
and (c) correspond to the system with hm = (phe- 
nomenological) and no lateral offset (k-p). Panels (b) and 
(d) correspond to the system with hm = 0.02 meV and 
3 nm offset, respectively. As highlighted by the dashed 
circles, the inclusion of finite hm and finite offset intro- 
duces anticrossings at the same positions. This strongly 
supports our conclusion that lateral offset is responsible 
for the phenomenological spin mixing term. 

A detailed derivation of the spin mixing term hm in 
terms of the k-p matrix elements induced by the misalig- 
ment potential can be found in Appendix [Bl It follows 
from the analysis that the proximity of LH states plays a 
critical role in the spin mixing of HH. The direct coupling 
between |3/2,1) and | — 3/2,2) states is small, because 
these two states are essentially HHs localized in opposite 
dots at the electric fields where the intersections occur. 
The coupling is mostly mediated by the excited states 
with Fz = ±1/2, which contain sizable LH components 
and are hence delocalized over the QDM for the entire 
range of electric fields under study. This allows them 
to couple the two "HH" states efficiently. A diagram 
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FIG. 5: (Color online). Exciton energy levels (dotted lines) 
and photoluminescence spectrum at B = 6 T calculated with 
Eq. ([2} (a-b), and k-p theory (c-d). Note the correspondence 
between (a) and (c), and between (b) and (d). The presence 
of QDM misalignment originates the same anticrossings as 
the spin mixing term km (highlighted by dashed circles). 



summarizing the coupling is shown in Fig. [51 To illus- 
trate the localization of the states, in-plane and vertical 
parts of the spinor envelope components arc written sep- 
arately. B (r) indicates localization in the bottom (top) 
dot, while {B ± T) indicates bonding and antibonding 
delocalized states. Mixing occurs between states with 
the same Jz- The coupling between the "HH" and delo- 
calized "LH" states is moderately strong. The hole spin 
mixing is limited by the coupling between the two delo- 
calized "LH" states because these two spinors accumulate 
weight on opposite components. 

To close this section we mention that the existence of a 
lateral offset between the QDs forming the QDM is quite 
plausible. In the right panel of Fig. [3l we show cross- 
sectional scanning tunneling microscopy (XSTM) data of 
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FIG. 6: (Color online) Diagrammatic representation of the coupling between two localized "HH" states with opposite pseudo- 
spin and location, as indicated by the single particle notation. All four spinor components of each \F^,k) quantum state are 
shown, with the shaded boxes indicating the approximate relative weights of each component. The interaction between localized 
'HH' states is mediated by delocalized states with large LH components. Note that the spinor labels are only approximate: 
misalignment breaks the symmetry, mixing different spinor configurations and enabling hole spin mixing. 



vertically stacked InAs QDs showing a lateral offset along 
the stacking axis. It is difficult to get good statistics on 
the occurrence of lateral offsets from XSTM data because 
it is sensitive only to lateral offsets parallel to the cleavage 
plane. A survey of our XSTM data indicates that lateral 
offsets are not uncommon in our QDMs. However, the 
lateral offsets are typically much smaller than the 5 nm 
offset found to fit the spectra presented in Fig. O This 
is consistent with the observation that the degree of hole 
spin mixing in this example is unusually large. 



V. SUMMARY 

We observe optical intensity from dark excitons in a 
number of QDM samples. This dark exciton optical in- 
tensity can be phenomenologically explained by the pres- 
ence of hole spin mixing. We have presented an exam- 
ple in which the hole spin mixing is sufficiently large 
and the tunneling sufficiently small that bright-dark exci- 
ton anticrossings can be directly observed and measured. 
We have used a phenomenological matrix Hamiltonian to 
demonstrate that the inclusion of hole spin mixing qual- 
itatively explains all of the new anticrossings in the ex- 
perimental data. We have then used a k-p theory to show 
that molecular symmetry breaking in the form of QD lat- 
eral offset changes the mixing of HH and LH states and 
enables the same form of hole spin mixing. The combina- 



tion of experimental evidence and k-p theory lead us to 
the conclusion that hole spin mixing can occur in QDMs 
as a result of symmetry breaking. 

A quantitative understanding of the spin mixing term 
will require detailed experiments in which the molecu- 
lar symmetry can be broken in a quantifiable way. It 
will also require the development of more sophisticated 
theory and modeling techniques to quantitatively deter- 
mine the magnitude of spin mixing and its dependence 
on molecular structure and symmetry. The experimental 
evidence and k-p theory we present suggest that there is 
a rich regime of spin physics to be explored in QDMs. In 
the final section we discuss the implications of hole spin 
mixing for device applications. 



VI. IMPLICATIONS OF HOLE SPIN MIXING 

As discussed earlier, the low-lying hole states with 
~ ±3/2 have a strongly dominant HH component 
and hence can be assigned a pseudo-spin. However, at 
the electric fields where the anticrossings induced by mis- 
alignment occur, these states mix with other spinors with 
opposite pseudo-spin. To quantify the effect of such mix- 
ing on the spin purity, in Fig. [7| we plot the expectation 
value of the hole pseudo-spin Sh for the first and second 
molecular excited states of Fig. U] (i.e. the states that 
arc a mixture of = +3/2,1) and -3/2,2)) 
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as a function of the applied electric field. {Sh) is evalu- 
ated from the weight of the components of the Luttinger 
spinor as 



{Sk 



/3 - 



(5) 



where the factor 1/3 acting upon the LH components 
comes from the Bloch function coefficients. Fig. [7^ and 
Fig. [Tb correspond to the QDM without and with lat- 
eral offset, respectively (as in Figs. 2^ andlDs). One can 
see that in the absence of misalignment (Sh) is nearly 
pure (±1/2)4^ but the inclusion of misalignment severely 
degrades the spin purity at electric fields near the spin 
anticrossing points. 



1/2 



00 
V 



-1/2 
1/2 



M 

V 



-1/2 



■ 1st excited state 

■ 2nd excited state 



-4 -2 2 4 

F (kV/cm) 

FIG. 7: (Color Online) Expectation value of the hole pseudo- 
spin as a function of applied electric field. Red and blue lines 
correspond to the first and second excited states of Fig. [l] (a) 
no lateral offset, (b) 3 nm lateral offset. The reversal of red 
and blue lines in panel (a) arises when the energy order of 
the two states reverses and does not indicate a degradation 
of hole spin purity. 

These perturbations in the hole spin purity further il- 
lustrate that the HH-only approximation has only limited 
validity in QDM. Near the anticrossing points, the effect 
of valence band mixing must be taken into account when 
designing information storage or manipulations protocols 
based on the hole spin projections. The contributions of 
the other spinor components may enable additional de- 
coherence or dephasing mechanisms that limit the ability 
to store and manipulate quantum information, as well as 
the preparation of pure spin hole states through the pos- 
itive trion^i^. These mechanisms could be suppressed 
by minimizing the structural distortions that break the 
molecular symmetries or by designing spin storage and 
manipulation protocols that are insensitive to the mix- 
ing of additional hole spin components. Away from the 



anticrossing points, the HH spin projections remain rela- 
tively pure and viable for quantum information storage. 

The spin mixing mechanism we describe here is rem- 
iniscent of that reported by Ferreira and Bastard for 
single asymmetric quantum wells Holes in symmetric 
quantum wells have a well defined pseudo-spin in spite of 
the HH-LH coupling because the parity symmetry pre- 
vents undesired mixing between spinorsi^^^ Introduc- 
ing envelope function asymmetries that break the parity 
bring about a D'yakanov-Perel-like mechanism of spin 
relaxation]^ In our QDM system, the vertical parity 
symmetry is lifted by the lens-shaped confinement and 
the different composition of the dots. Still, this does 
not suffice to mix pseudo-spins because the lateral con- 
finement, which is obviously absent in quantum wells, 
imposes an additional symmetry, namely circular sym- 
metry. As a result, the "spin up" and "spin down" hole 
states have different total angular momenta and they 
remain orthogonal. In order to mix the orthogonal spin 
projections, one has to break the rotational symmetry. 
As we have shown, this is most efficiently achieved by 
the misalignment of the QDM. 

The envelope origin of the SO term we report sug- 
gests that it can be controlled with external field, and 
hence used in a similar fashion to the Rashba term for 
spintronic applications Indeed, the magnitude of the 
spin anticrossing gaps observed in Fig. [Ife) are compa- 
rable to those of electrons in InAs QDMs with strong 
Rashba interactioui^ This makes holes in QDMs par- 
ticularly suitable for electric-field-induced spin manipu- 
lation. For example, two orthogonal spin states of the 
hole in the bottom QD could provide the qubit basis. 
Away from the electric field of tunnel coupling and hole 
spin mixing, these states are well isolated and could pro- 
vide a robust means of storing quantum informatioui^ 
The applied electric field could then be varied to bring 
the dots into the regime where hole spin mixing becomes 
strong in order to mix the two spin configurations. A de- 
tailed examination of the interaction strengths and elec- 
tric field pulse sequences required to effect a single qubit 
rotation are beyond the scope of this paper. If feasible, 
however, such a scheme would require fewer resources 
to implement spin control protocols than current pro- 
posals. Spin rotations could be implemented with only 
a single applied electric field, eliminating the need for 
pulsed lasers tuned to transitions specific to each dot or 
GHz frequency fields used to implement g-TMR rotations 
of spin projectionsi^i^'^ At the same time, the QDM 
structure preserves the opportunity to use optics for spin 
initialization, readout and control of 2-qubit operations. 
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APPENDIX A: EMPIRICAL PARAMETERS 

The numerical values of the parameters used in the 
matrix Hamiltonian are not free fitting parameters. The 
value of each parameter can be determined from experi- 
mental data, txo = —0.107 meV is determined from the 
measured zero magnetic field anticrossing. g^s = —1-695 
is determined from the Zeeman splitting of the lower en- 
ergy line well away from the anticrossing region. Simi- 
larly, ghT = —1.66 is determined from the Zeeman split- 
ting of the higher energy indirect line well away from 
anticrossing. ge = —0.6 and 6o = 0.101 meV are deter- 
mined from the asymptotic energies of the dark states. 
d = 5.8 nm is determined from the slope of the indirect 
transition energy. B = 6 T is the known value of the 
applied magnetic field, /is = 0.0579 meV is the Bohr 
magneton. 1712 = 0.47 is fit to the data by looking at the 
higher energy anticrossing, which shows little affect from 
the hole mixing, hm ~ 0.092 meV in the calculations of 
FiglJt and f is fit by looking at the new bright-dark anti- 
crossings. The suppression of txo in the calculations with 
nonzero hm is intuitively expected because the lateral off- 
set increases the distance between the center of the QDs. 
k-p calculations confirm that the LH influence does not 
alter the suppression of txo ■ Note again that the calcu- 
lations in Figl5]use both electron spin projections, so the 
matrix and intensity vector analogous to Eqn.Q for the 
electron spin up case is also used. A value of Ind = 0.5 
is used to generate the calculated spectral maps. 



APPENDIX B: THEORETICAL MODEL AND 
K P DERIVATION OF SPIN MIXING TERMS 

In this Appendix we describe the theoretical model 
employed for the k-p calculations and then derive an ex- 
pression for the spin mixing operator of Eq. ([2]), hm, in 
terms of the matrix elements induced by misalignment. 

We write the Hamiltonian of the QDM as: 

= WlK + VqdM + Kffsct • (Bl) 

Here TYlk is the three-dimensional Luttinger-Kohn 
Hamiltonian, including longitudinal magnetic and elec- 
tric fields. Vqdm is the confinement potential of an ideal 
QDM, formed by two vertically stacked lens-shaped QDs 
(spherical caskets) with perfect alignment (see Fig. |5^) . 
The potential is zero inside the dots and Vo outside. 
K)ffset is a perturbative potential induced by laterally off- 
setting the QDs. It is the difference between the potential 
of misaligned and ahgned QDM (see Fig.IBb). 

We first obtain the eigenstates of Hsym ~ T^-lk + Vqdm 
and then project the full Hamiltonian Ti. into these states. 



a b 




FIG. 8: (Color Online) (a) Confinement potential of an ideal 
QDM. (b) Perturbative potential introduced by QDM mis- 
alignment. Blue and red regions correspond to +Vo and — Vo, 
respectively. 



The minimal basis set which captures the spin mixing fea- 
tures of Eq. ^ is formed by the six lowest-lying eigen- 
states of Tisym, namely: 
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The above states are obtained by numerically integrat- 
ing the Luttinger-Kohn Hamiltonian, written in cylindri- 
cal coordinates, with a finite difference scheme.— The 
complete Hamiltonian Ti. is then integrated using an ex- 
act diagonalization technique. The QDM is constructed 
of two QDs with radius 15 nm, height 2 nm and an inter- 
dot separation of d = 1.7 nm (unless otherwise stated). 
We use InAs Luttinger parameters 71 = 20.0, 72 = 8.5 
and 73 = 9.2i^ The valence band offset is Vo = 200 meV, 
and a constant g-factor gh = —1.5 is assumed. 

The energy of the Fz = ±3/2 states as a function of the 
electric field is represented in Fig. 2^ for an ideal QDM at 
B = 6 T. The F^ = ±1/2 states are a few meV higher in 
energy. From the dominant component of the spinor, one 
can identify the states |1) and |2) (|5) and |6)) with the 
-IJ. (ff) HHs of Eq. 11]). These states are localized cither in 
the top or the bottom dot, except for a narrow window 
near the resonant field Ez ^ where they form delocal- 
ized bonding and antibonding states. On the other hand, 
the excited states |3) and |4) have a significant admix- 
ture of HH and LH components. This allows them to 
be significantly delocalized for all the values of Ez in the 
figure, forming bonding molecular states4i 

The exciton calculations of section ITVl are carried out 
using the hole states calculated as described above and 
electron states calculated in a similar fashion but with a 
single-band effective mass model<^ The electron mass is 
771* = 0.06, the conduction band offset is set to Vo = 500 
meV and the g-factor = —0.6. To mimic the exper- 
imental situation, we force the electron to stay in the 
bottom QD (a single dot potential is used). Electron-hole 
Coulomb interaction is accounted for using a configu- 
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ration interaction method on the basis of the Hartree 
products formed by the electron ground state and the six 
lowest hole states. The matrix elements are integrated 
using Monte-Carlo routines, and the exciton emission 
intensity is computed within the dipole approximation!^ 

We can obtain a pcrturbativc expression for the spin 
mixing term hm by projecting Hamiltonian (|Bip into the 
basis of the states |1) to |6). This yields: 



This can be understood from an analysis of the matrix 
elements. For example, (l|V^ffsot|5) reads: 



<3/-3(r)(J. = +|| c+3/o(r)|J,^+|) 

c\ /_2(r)(J. - +il „ c_ . A(r)| J. ^ +1) \ 
<:/-i(r)(J. = -i| c,./.(r)|J. = -i) /' 

c' 3 /+o(r)(J. = -II c_3 h{r)\.h - -I) 



( El Vi2 V^i3 \ 

^12 E2 1/23 
^1*3 ^2*3 ^34 

T/45 Vie 

^4*5 E, ^56 

V Vie n*6 Ee J 

where Ei is the energy of the hole state while Vij = 
(* I Vl^ffset |j) is the matrix element induced by the misalign- 
ment potential, which couples the states \i) and 

For simplicity, in Eq. (jB2[) we have set to zero the 
matrix elements which are small or non relevant to the 
spin mixing under discussion. In particular, we note that 
there is no direct mixing between the HH JJ. and states. 



where I is the identity matrix. While |1) gathers most 
of its weight in the Jz — —3/2 component, which has the 
lowest envelope angular momentum (ruz = 0), |5) does so 
in the Jz = -1-3/2 component. Since each component of 
the bra couples to that of the ket with equal Jz , the cou- 
pling between the two vectors is negligible. By contrast, 
both |1) and |5) couple with the excited states |3) and |4). 
This is because such states have a strong admixture of 
components and arc delocalized over the entire QDM, so 
they are able to couple with "HH" states of cither "spin" 
regardless of their localization. 

In order to compare with Eq. we reduce Hamilto- 
nian (|B2|) to an effective 4x4 Hamiltonian on the basis of 
"HH" states only (i.e. states |1), |2), |5) and |6)). After 
some algebra, we obtain: 



/ Ei + k{X - Ei)Vi3V{:, V12 + fc(A - £;4)1^13^2*3 

^1*2 + fc(A - £^4)^1*3^23 E2 + fc(A - Ei)V23Vi:i 

kV^AV^^Vir, kVuV^%V45 E, + kiX^E3)V45V:, l^e + A:(A - £3)145^6 

\ kV34V*3Vie fcl/34 1^2*3 ^46 ^56 + k{X - £^3)^4*5 ^46 Ee + k{X - ^3)146^*6 / 



fc^34^13^4 



45 



^^34^23^5 



^^34^13^6 
^^34^23^6 



fB3) 



where A arc Hamiltonian (|B2p eigenvalues and k = 

iiX-E3)iX~Ei)-V3iViir\ 

Hamiltonian (jBSp can be simplified if we disregard the 
imaginary part of the matrix elements, which is neglegi- 
ble in our QDM because K>ffset has even parity along the 



azimuthal direction, and we further consider that the or- 
bital part of |1) and |5) (|2) and |6)) are very similar, 
as they differ in the pseudo-spin only. It follows that 
V12 « V56 = a, Vi3 « V45 = a and ^23 ~ ^45 = /3. The 
Hamiltonian can then be rewritten as: 



/ Ei + k{X - Ei)a 
a + k{X - Ei)al3 
hmi 

\ hm 



2 a + k{X- E4)aP 
E2 + k{X - Ei)l3'' 



hm 
hm2 



hmi 
hm 

E5 + k{X - E3)a'^ 
k{X - E3)af3 



hm \ 

hm2 

a + k{X- E3)af3 
Ee + k{X - E3)P^ J 



(B4) 



where hm = kV^ia/], hmi = kVs^a^ and hm2 = kV^iP^. 

By inspecting the effective Hamiltonian (jB4[) . a num- 
ber of conclusions can be drawn. First and foremost, wc 



have found an expression for the spin-mixing term hm 
in terms of misalignment matrix elements. This clearly 
reveals the importance of the indirect coupling through 
the excited \Fz = ±1/2) states (V34). Indeed, hm is in- 
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versely proportional to the energy splitting between the 
\F^ = ±3/2) and the excited \F^ = ±1/2) states (through 
k). Second, in addition to hm there are other spin mix- 
ing terms, namely hmi and hm2. These terms couple 
"HH" states localized in the same dot. Fig. IB] shows an 
schematic representation of the low-energy hole states 
of a QDM at finite magnetic field and the spin mixing 
terms coupling each pair of states, hmi and hm2 are 
not needed in Eq. ([2]) because the states they mix are 
split by the Zeeman energy, which in the B = 6 T ex- 
periment is far larger than the SO perturbation. Finally, 
we note that misalignment not only induces spin mixing, 
but also affects the "HH" energy levels (through the di- 
agonal terms of (jB4p ) and the spin-conserving tunneling 
rates. This can be observed in the k-p results in Fig. H) 




Electric field 



FIG. 9: Schematic representation of the hole energy levels of 
Fig. |4] indicating the level mixing induced by the spin-orbit 
terms hm, hmi and hm2 of Hamiltonian (|B4|) . Solid (dashed) 
lines are used for states with pseudo-spin JJ. (-ff). 
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